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Abstract 

The phonon damping model (PDM) is extended to include the effect of angular momentum at 
finite temperature. The model is applied to the study of damping of giant dipole resonance (GDR) 
in hot and noncollectively rotating spherical nuclei. The numerical results obtained for 88 Mo and 
106 Sn show that the GDR width increases with both temperature T and angular momentum M. 
At T > 4 MeV and M < 60 h the increase in the GDR width slows down for 106 Sn, whereas at 
M > 80h the GDR widths in both nuclei nearly saturate. By adopting the nuclear shear viscosity 
extracted from fission data at T = 0, it is shown that the maximal value of the angular momentum 
for 88 Mo and 106 Sn should be around 46 and 55 h, respectively, so that the universal conjecture 
for the lower bound of the specific shear viscosity for all fluids is not violated up to T = 5 MeV. 
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I. INTRODUCTION 



The giant electric dipole resonance (GDR) in highly excited nuclei has been a subject of 
a considerable number of theoretical and experimental studies over the last three decades. 
At present, comprehensive experimental data on the dependence of the GDR width on 
temperature T and angular momentum J have been accumulated for a large number of 
medium and heavy-mass compound nuclei starting from potassium up to lead isotopes (See 
Ref. [l] for the recent compilation). This systematic shows that the full width at half 
maximum (FWHM) of the GDR increases significantly with T at low and moderate T, but 
seems to saturate at high T, above T ~ 4 - 5 MeV. At a given value of T, the GDR width 
also increases with J, but this increase is noticeable only starting from J > 25 h for copper 
and heavier isotopes. 

A number of theoretical approaches has been proposed to describe the behavior of the 
GDR width as a function of T and J and two monographies have been published on this 
subject {2, 3|. Among the theoretical approaches to the study of the hot GDR, the phonon 
damping model (PDM) {4- 7] offers a microscopic mechanism of the GDR damping via 
coupling of the GDR to non-collective particle-hole (ph) transitions, which exist already at 
T = 0, and particle-particle (pp) ones, which appear at T ^ because of the distortion 
of the Fermi surface. The PDM is able to describe consistently both the increase in the 
GDR width at low and moderate T as well as its saturation at high T by using the single- 
particle energies, obtained in the Woods-Saxon potentials for protons and neutrons, and a 
set of 3 parameters, fixed at T = 0, which are the GDR energy before the coupling, the 
ph and pp coupling constants (4 -6]. By including nonvanishing thermal superfluid pairing, 
the PDM is also able to explain a nearly constant value of the GDR width at T < 1 MeV, 
found in the inelastic scattering of 17 on 120 Sn In a recent development, the PDM was 
employed to calculate the ratio rj/s of the shear viscosity rj to the entropy density s directly 
from the GDR photoabsorption cross sections in medium and heavy spherical nuclei at T ^ 
js]. The results obtained show that this ratio, which is called the specific shear viscosity 
hereafter, decreases with increasing T to reach (1.3 — 4) KSS at T = 5 MeV. The quantity 
KSS = h/(4irkB) has been conjectured to be the lower bound for all fluids (9)]. 

A shortcoming of the PDM is that it does not include so far the dependence on angular 
momentum. The aim of the present paper is to remove this deficiency by extending the PDM 
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to incorporate the effect of angular momentum to make the model capable to describe the 
dependence of the GDR width on both T and J. The present extension will also be applied 
to examine the credibility of the preliminary analysis of the recent data [loj], according to 
which the GDR width extracted at T ~ 4 MeV and J = 44 h seems to be smaller than that 
measured at T ~ 3 MeV and J = 41 h in 88 Mo. 

The paper is organized as follows. The PDM is extended to include the effect of angular 
momentum in Sec. [Ill The analysis of numerical results is presented in Sec. IIII1 The article 
is summarized in the last section, where conclusions are drawn. 



II. ANGULAR-MOMENTUM EFFECT WITHIN THE PDM 
A. Extention of PDM to finite angular momentum 

The axis of quantization (or the laboratory-frame z axis) of a spherical nucleus can always 
be chosen to coincide with the body-fixed axis. The latter is aligned with the direction of the 
total angular momentum within the quantum mechanical uncertainty 
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-Il3|. As the result, 



the total angular momentum is completely determined by its ^-projection M alone. The 
rotation of a spherical nucleus (as of a classical system) about the z axis is called noncollective 
or single-particle rotation. On the other hand, even for non-spherical nuclei, especially the 
axially symmetric ones, at high excitation energies and/or high T, the values of angular 
momentum projection on the symmetry axis are mixed among the levels in the region of 
high level densities, which worsen the axial symmetry. The melting of shell structure will 
also eventually drive nuclei to their average spherical shape. The PDM is extended here to 
include the angular momentum effect for such noncolleticve rotation at finite T. 

For this purpose, the Hamiltonian of a spherical system, which is rotating about the 
symmetry z-axis, can be written in the following form 

H = H - 1 M J (1) 

where H is the PDM Hamiltonian described in Ref. , and M represents the total 

angular momentum J, which, in this case, coincides with its z-projection M. The latter is 
defined as 

M = J2 m *( N k-N- k ) , (2) 

k>0 



where, for simplicity, the subscripts k are used to denote the single-particle states \k,mk) 
in the deformed basis with the angular momentum k and the positive single-particle spin 
projection m k , whereas the subscripts —k denote the time-reversal states \k, — mj~) (m^ > 
0). In the spherical basis this corresponds to \j,m). The deformed basis is preferable here 
because the rotation about the z- axis eventually resolves the (2j + 1) degeneracy of each 
spherical orbital j. The particle number operator N consists of 

N = ^(iV fe + N. k ) , N ±k = a ] ±k a± k , (3) 

fc>0 

where a* ±k (a± k ) denotes the creation (annihilation) operator of a particle with angular 
momentum k, spin projection ±mfc, and energy e k . 

By using Eqs. §2§ and ([3]), the Hamiltonian transforms in the deformed basis into 

H = ^(e k - A - ^m k )N k + ^(e fc - A + ^m k )N_ k + ^ u qQ\Qq 

k>Q k>0 q 

+ E + aV-vXQj + Qg) , (4) 

k,k'>0 q 

where A denotes the chemical potential. The particle (p) states are defined as those with 
e k > A, whereas the hole (h) states are those with e k < A. The operator Q^ q (Q q ) creates 
(annihilates) a phonon with energy u q , which describes the collective vibration. Therefore, 
the right-hand side of Eq. (j2J) describes two mean fields, one of the single particles as the 
first two terms, another of the phonon field associated with the GDR as the third term, and 
the coupling between them in the last term with matrix elements J 7 ^,- This coupling causes 
the damping of the GDR. The effects of angular momentum, which are incorporated in the 
first two terms, resolve the degeneracy of spherical orbitals. As the result, each of spherical 
orbital j with energy ej splits into 2Qj — 2j + 1 distinctive levels, half of which consists 
of levels with energies e k + 'yrrik , whereas the other half consists of levels with energies 
e k — 7wifc, with k — 1, £1/2, where £1 = 2 Y^j £lj is the total number of levels. Thermal 
superfluid pairing plays the role of keeping the GDR width essentially unchanged at T < 1 
MeV at M = [7]. Pairing has a negligible (or no) effect on the GDR width in the region 
of moderate (high) T and M, which is the main focus of the present paper. Therefore it is 
not included in the calculation of the GDR width here for simplicity. 

The chemical potential A and the rotation frequency 7 are two Lagrangian multipliers, 



which are defined from the equations for conservation of angular momentum M 



rjJ4l3| and 



particle number N as 

M = J2 m ^ ~ fk) , N = J2(fk + fk) . ( 5 ) 

k k 

where M = (M), N = (N), f± = (N± k ) with (...) = Tr[...exp(-/3if)]/Tr[exp(-/3#)] 
denoting the grand canonical ensemble average with /3 = T _1 . The PDM calculations at 
M = showed that, although the particle-phonon coupling in the last term of Eq. (jlj) is 
sufficient to cause the GDR width of around 4 to 5 MeV for medium and heavy nuclei at 
T = 0, it is not sufficiently large to cause a noticeable damping of each single-particle state. 
Therefore, the single-particle occupation numbers f k can be well approximated with the 
Fermi-Dirac distribution for nonineracting fermions, namely 

= exp(^) + 1 ' ^ = c *- A= F7m fc . (6) 

By solving Eqs. (jSJ) at each value of temperature T, the chemical potential A and rotation 
frequency 7 are defined as functions of T. Following the same procedure as that employed 
to obtain the PDM equations at M — 0, one now proceeds to derive the strength function 
and the width of the GDR atM^O below. 



The following double-time retarded Green functions are introduced, which describe 
1) The propagation of the free phonon: 



G q (t-t') = ((Q q (t);Ql(t'))) , (7) 

2) The transition between a nucleon pair with positive spin rrik and the phonon: 

Gtk' q (t-t')-((a{(t)a kl (ty,Ql(t'))) , (8) 

3) The transition between a nucleon pair with negative spin —m k and the phonon: 

W* - = ((aUt)a-k>(t); QW))) , (9) 

where the standard notation G(t,t') = ((A(t); B{t'))) = -%9{t - t')([A(t), B{t')}±) is used 
for the double-time retarded Green function with [A(t),B(t')]± = A(t)B(t') - B(t')A(t) for 



boson operators, and A(t)B(t') + B(t')A(t) for fermion ones 14j. Applying the standard 
procedure of deriving the equation of motion for the double-time Green function in the form 
idG/dt = 8(t -t')([A(t),B(t)] ± ) + (([A(t),H(t)];B(t'))) to the Green functions © - © 
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and Hamiltonian (j3J), one obtains an hierarchy of equations. This infinite series contains, 
beside the Green functions ([7]) - (jUJ), also the higher-order ones. To close this set, a simple 
decoupling is used by pairing off, wherever possible, the operators referring to the same 
time, e.g. {(a± k (t)a±k'(t)Q q (t)] Q\t'))) = 5kk'ftG q {t — t'). The whole procedure is the same 
as has been described thoroughly in Refs. |4j, [5| so it is not repeated here. The final set 
contains only 3 coupled equations for three double-time Green functions ((Tj) - Q, whose 
Fourier transforms into the energy variable E are 



(E - u q )G q (E) = -L + J2 4 q Mk> q (E) + Gw q (E)\ 

kk' 



(E - E- + E-)g+ lq (E) = - f+)G q (E) , 

(E - Et + E+)g- k ,(E) = - f k )G q (E) . 



(10) 

(11) 
(12) 



Expressing G kk r q (E) in terms of G q (E) by using the last two equations, and inserting the 
results into Eq. ( TLOl) . one obtain the final equation for the Green function G q (E), which 
describes the phonon propagation, as 



G q {E) 



2n E — uj n 



+P q (E) , P q (E) = 



(9)12 



kk' 



fk' fk 



fk' fk 



E ~ E k + E k , E — E k + E k , 

(13) 

The principal value of the polarization operator P q {ui) at a real u defines the energy shift 
from the unperturbed phonon energy u q to u q under the effect of particle-phonon coupling 
in the last term of the Hamiltoanian (j3J), whereas the phonon damping 7 g (cu) is defined 
as the imaginary part of the analytic continuation of P q {E) into the complex energy plan 
E = co ± ie, that is %(u)) = QmP q (uj ± ie). The final results reads 
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T,F&n(ft - tiM" - E k + Ey) + Uk' - KMu - Et + 25+)] , (14) 



kk' 



which, by using the ^-function representation 8(x) = \im e ^ e/[n(x 2 + e 2 )}, transforms into 

fk'~f'k 1 fk' ~ fk 



kk' 



[CO 



E^ + E ki y + e 2 



+ 



[u - Et 



E+y + e 2 



(15) 



J k ^ ^k'J ^ c ^k ^ ^k> 

The spectral intensity J q {oS) is found by definition from the analytic properties of the 
Green function G q (E) (TI3"]) as 

.G q (u + ie) - G q (oo -ie) 1 j q (u)(e^ - 



J q {u) 



7T (U - UJ q ) 2 + tf(u) ' 



(16) 



The GDR strength function S(u>) is obtained from the spectral intensity J g (u>) as S(co) = 
J q (u)[exp((3uj) — 1], where J q (uj) denotes J q (to) calculated at the GDR energy oj = Egdr(T) 
at temperature T [#46| • The result reads 

OA A = I llM (1 7 ) 

The FWHM width T(T) of the GDR is defined as a function of T at each value M of the 

fin 

total angular momentum as |4j, [5] 

T(T) = 2 lq [u = E GDR (T)} . (18) 

This width consists of the quantal and thermal parts [See Eqs. (la) - (lc) of Ref. Q]. The 
quantal part is the spreading width T^, caused by coupling of the GDR to ph configurations. 
The thermal part comes from coupling of the GDR to pp and hh configurations, which 
appear due to the distortion of the Fermi surface at T / 0. The escape width T^, which 
arises due to coupling to continuum and is related to the direct decay by particle emission, 
is usually small (in the order of few hundreds keVs) [15], and there is no evidence that it 
is sensitive to the change of T in medium and heavy nuclei. In the numerical calculations 
within the PDM, the effect of the escape width is taken into account via the smoothing 
parameter e in Eq. (EES}, which usually does not exceed 1 MeV (See Sec. IIHIB). 

Finally, it is worth mentioning that Ref. 16) proposed an evaporation width T ev of 



the compound nucleus states to be added twice to the total GDR width (See also p. 226 
of Ref. [2 1). This evaporation width comes from the quantum mechanical uncertainty 
principle, according to which, the energies of the compound nucleus with a finite lifetime t 
cannot be known with an accuracy better than T ev ~ H/t. Once this width is incorporated 
in the calculations of the transition probabilities that define the amplitudes of the GDR 
photoabsorption cross section, it gives the natural width of 2r e „, which is approximately 
the sum of the widths of the initial and final states. In heavy nuclei at high excitation 
energies, the width 2r e ,; can amount to few MeVs (See, e.g., Tabs. 3 and 4 as well as Figs. 
5 and 8 of Ref. jl3])- The PDM does not takes the evaporation width T ev into account. 
However, it should be noticed that, a complete procedure requires the comparison of not 
just the GDR energy and width, but of the entire GDR shape generated by the theoretical 
strength functions with that obtained from the measured 7-ray spectra. For this purpose, 



the authors of Ref. 



18) proposed a method, which incorporates the theoretical strength 



functions directly into all the decay steps of the full statistical calculations. This method 
allows to test the contribution of the evaporation width as well. As a matter of fact, by using 
the evaporation width 2T ev obtained in Refs. [if], [l7|, the authors of this method found that 
the contribution to the total 7-spectrum by the evaporation with is small relatively to the 
total spectrum including all the decay steps. The overall high energy 7-ray spectra resulting 
from the complete CASCADE calculations including the evaporation width 2T ev turn out 
to be essentially identical to those obtained by neglecting this evaporation width even up 
;o the excitation energy higher than 120 MeV for 120 Sn (i.e. at T > 3.3 MeV). In Ref. 



191 ] . by using the same method, the PDM strength functions were also incorporated into all 
the steps of the CASCADE calculations, and a reasonable agreement between theory and 
experiment were found, especially at T > 2 MeV. These results indicate that the effect of the 
evaporation width 2T ev may become substantial only at much higher values of temperature 
and angular momentum (T ^> 3.3 MeV and J ^> 30h). Further experimenal tests of this 
effect at high excitation energies are still required. 

B. Specific shear viscosity 

Knowing the GDR parameters Ecdr{T) and T(T), one can calculate the GDR photoab- 
sorption cross-section <jgdr(u,T), and the shear viscosity 77 (T) as functions of T at each M 



by using t 
(See Ref. 



re Green-Kubo formula in combination with the fluctuation-dissipation theorem 



8] for details). The result yields [See Eq. (6) of Ref. jg]] 

= lim W.r) = ( ) rm WQ)' + [r(0)/2]- (19) 

'> [1 > S C " {U> r(0) E GDR [Tf + [T(T)/2]= ' (1J) 

where the normalization parameter C is chosen to reproduce the value 77(0) of rj(T) at T = 
0, that is C = \\m^ [a GDR {u,T = 0)]/r/(0). 

To obtain the specific shear viscosity rj/s one needs to know the entropy density (entropy 
per volume V) s = S/V = pS/A with the nuclear density p = 0.16 fm -3 and nuclear 
mass number A. The entropy S at temperature T is obtained by integrating the Clausius 
definition of entropy as S = dr[T~ 1 d£/dT], where £ is the total energy of the system at 
temperature r. By using this definition and taking the grand-canonical-ensemble average of 
the PDM Hamiltonian (j3J), it follows that S = S a + Sq, where S a and Sq are the entropies 
of the quasiparticle and phonon mean fields, respectively. The quasiparticle entropy S a is 
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given in units of Boltzmann constant fcg as 

«5« = - 5>£W + ( X " " < ) + n k lnn k + (! " % M 1 " % )] ( 2 °) 

where = {exp TT^fc) + I} -1 are the quasiparticle occupation numbers on the fc-th 
level with quasiparticle energies Ef, = y/ (e^ — A) 2 + A(T). The pairing gap A(T) is needed 
in the calculation of the entropy in open-shell nuclei to ensure that S a goes smoothly to zero 
at T — Y 0. At T > T c ~ 0.57A(0), where the gap vanishes as in the BCS theory (or remains 
finite but small as in the case when thermal fluctuations of the pairing field are included), rvj: 
coincide with (or approach to) f k and Eq. ( 120]) yields the single-particle entropy. Because 
the effect of thermal fluctuations on the entropy is not large at moderate and high T, the 
present paper does not take them into account, adopting the pairing gap A(T) obtained 



from the BCS theory including the angular momentum effect at T ^ 111, Il2 | 



The phonon entropy Sq is given as Sq = + v q ) ln(l + v q ) — v q \nn q ], where v q ~ 

[exp(Pu q ) — l] -1 is the phonon occupation number approximated with the Bose-Einstein 
distribution for noninteracting bosons. However, in practical calculations of GDR within 



the PDM, Uq is close to Egdr(T) 3> T, which means that v q are small [8(. As the result, Sq 
for GDR is negligible as compared to S a , allowing Sq to be safely omitted to give S ~ S a 
in the present calculations. 

III. ANALYSIS OF NUMERICAL RESULTS 
A. Ingredients of numerical calculations 

The single-particle energies €j are obtained from the spherical Woods-Saxon potentials for 
neutrons and protons. At M — 0, these (2j + l)-degenerate levels span a large space starting 
from the bottom lsi/2 level located at around —40 MeV for neutrons and —30 MeV for 
protons up to around 22 MeV for 88 Mo and around 26 MeV (18 MeV) for protons (neutrons) 
in 106 Sn. They are kept unchanged with T based on the estimation within the temperature- 
dependent self-consistent Hartree-Fock calculations 20], which have demonstrated that the 



single-particle energies in heavy nuclei weakly change with T up to T ~ 5 MeV. 

The PDM assumes that the matrix elements J 7 ^ for the coupling of the GDR to non- 
collective ph configurations, causing the quantal width already at T = 0, are all equal to 
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parameter T\, whereas those for the coupling of the GDR to pp (hh) configurations, Fpf and 
J 7 ^, causing the thermal width at T 7^ 0, are all equal to parameter JF 2 . The justification of 
such assumption has been discussed in detail previously so it is not repeated here [See, e.g., 
Sec. II B of Ref. [7j for the details]. With this assumption the PDM has three parameters, 
J- 1, J-2, and the unperturbed phonon energy u q . The latter, in the case of GDR (q = 1), is 
chosen to be close to the E GDR (0). The parameters T\ and J-" 2 are adjusted to reproduce 
the experimental value of the GDR width at (T = 0, M = 0), and so that the GDR energy 
Egdr(T) does not change appreciably with T. At M — 0, E GDR (T) is usually determined 
from the equation E G £> R {T) = u q + P q [EcDR{T)]. At M / 0, to avoid solving this equation 
at each value of M, it is assumed, based on experimental systematics, that Egdr does not 
depend on T and M. Therefore, the values of parameters T\ and Ti, chosen at M — 0, are 
adopted, whereas the experimental value of E GDR {T) = E GDR (0), which is 15 MeV for 88 Mo 
and 15.5 MeV for 120 Sn, is used in Eq. ( 117)) as T and M vary. 

For the calculations of the entropies in 88 Mo and 106 Sn, the pairing interaction parameters 
are chosen to reproduce the experimental values of the neutron and/or proton gaps at T = 
and M — 0. Hence, the present paper adopts Gz = 0.154 MeV for protons and Gn = 
0.153 MeV for neutrons in 88 Mo to obtain the pairing gaps A^(0) ~ 1.25 MeV and Atv(0) ~ 
1.39 MeV at M = 0. For the proton closed-shell 106 Sn the value G N = 0.124 MeV is used 
to give Ajv(0) ~ 1.25 MeV. These values of pairing gaps agree with the empirical data 21]. 
As for 77(0) in Eq. ( 1T9)) . the value 77(0) = lu, extracted from the GDR data at T = 22], 
and the lowest value 7/(0) = 0.6m, obtained by fitting the fission data with u = 10 -23 Mev s 
fm -3 are used [See Sec. II A of Ref. js| for the detail discussion on the selection of 7/(0)]. 



B. GDR Strength function and width as functions of T and M 

Shown in Figs. UJ and [2] are the strength functions S(u), obtained from Eq. (ITTj) . and 
the width, obtained from Eq. ( ITS)) , of the GDR in 88 Mo and 106 Sn at various values of 
temperature T and angular momentum M. A smoothing parameter e = 0.5 MeV in Eq. 
( TT5)) was used in calculations. The results of the GDR width (ITS)) do not change significantly 
within the interval 0.1 < e < 1 MeV, but the strength function S(u) is smoother by using 
a larger e. 

Three features are clearly seen from these figures. The first one is that the GDR width 
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FIG. 1: (Color online) GDR strength functions 5(w) for 88 Mo [(a) - (f)] and 106 Sn [(g) - (1)] at 
T = 0.5, 1, 2, 3, 4, and 5 MeV as shown at the curves in (a), and M = 0, 10, 20, 40, 60, and 80 H 
as shown in the panels. 
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FIG. 2: (Color online) FWHM of GDR for 88 Mo [(a), (c)] and 106 Sn [(b), (d)] as a function of T 
at several values of M (in h) shown at the curves [(a) and (b)], and as a function of M at several 
values of T (in MeV) shown at the curves [(c) and (d)]. The experimental data for GDR in 88 Mo 



. Q,y,u 



(triangles), 106 Sn (solid circles), 109 'H0g n (solid and open boxes) are adapted from Ref 
respectively. 

increases with both T and M. This increase is stronger at low T and M. The second feature 
is that, while the GDR shape becomes smoother as T and M increase, the smoothing caused 
by the angular momentum seems to be stronger than that caused by thermal effects (Fig. 
[Q. The third feature is that the GDR width approaches a saturation at moderate and high 
T or M. As a function of T, the width saturation begins at T > 4 MeV in 106 Sn [Fig. [2] (b)], 
whereas as a function of M the width saturation takes place in 106 Sn already at T > 3 MeV 
[Fig. [2] (d)]. In 88 Mo, which is lighter than 106 Sn, the GDR width keeps increasing with T 
up to T ~ 5 - 6 MeV at low M, but also starts to saturate at M > 60 MeV [Figs. M (a) 
and [2] (c)]. At very high M (above 100 h) the GDR width ceases to change with T, but this 
value of M is too large to be realistic as will be discussed in the next section. The width 
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FIG. 3: (Color online) (a) Rotation frequency 7 as a function of angular momentum M at several 
temperatures T = 0.4 (red solid), 1 (green dotted), 2 (purple dashed), and 4 (blue dot-dashed) 
MeV in 106 Sn. (b) Differences f k , — f k for the pair of neutron levels (W5/2 — 2>Pi/ 2 ), (W5/2 ~~ 2d 5 / 2 ), 
and (I59/2 — 2/7/2) in 106 Sn, calculated at T = 4 MeV, as functions of M. The signs (+) and (-) 
at the notations for the lines correspond to (f k , — f k ) and (f k , — /jT), respectively. 

saturation at high T is a feature, which is obtained within the mechanism of GDR coupling 
to pp (hh) configuration at T 7^ and has been discussed previously (See Refs. (4 , 5 ] , e.g.). 
The same mechanism takes place at M 7^ 0, where the contribution of GDR coupling to pp 
(hh) to the width appears whenever f k , — f k is not zero [(k, k') = (p,p') or (h, h')] [See Eq. 
( IT5|) ]. At small T the rotation frequency 7, which dictates the dependence the GDR width 
on M, clearly increases with M at low M but saturates at high M as shown in Fig. |3] (a) for 
106 Sn at T = 0.4 and 1 MeV. The temperature dependence of the GDR width is determined 
essentially by that of the differences f k , — f k at the right-hand side of Eqs. ( TT4l) and (|T5]l . 
An example of such dependence is shown in Fig. [3] (b), where the factors f k , — f k for the 
pair of neutron levels (ld 5 / 2 — 2pi/ 2 ), (ld 5 / 2 — 2d 5 / 2 ), and (lgg/ 2 — 2/7/2) in 106 Sn, calculated 
at T = 4 MeV, are plotted as functions of M. It is clear from this figure that f k , — at a 
given T approach a saturation at high M. This means that, although at large T the rotation 
frequency 7 increases with M even at high M [See Fig. [3] (a) at T > 2 MeV], the width 
saturation caused by temperature seems to dominate, which leads to the combined effect of 
width saturation at high T and M. 

Few experimental data, available only within narrow regions of T and M, are also shown 
in Fig. [2j These are the GDR widths in 88 Mo [10] and 106 Sn [24], as well as the neighboring 
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FIG. 4: (Color online) (a) - (c): Pairing gaps as functions of T at different values of M (in K) as 
indicated in the panels for protons (a) and neutrons (b) in 88 Mo and neutrons in 106 Sn (c). Panels 
(d) and (e) show the entropies as functions of T in 88 Mo and 106 Sn, respectively at several values 
of M. 



isotopes 109 ' 110 Sn 25j, which are also spherical in the ground state (i.e. at T =0 and M — 0). 
Except for the data point in 88 Mo at T = 4 MeV and M = 44 h [Figs. |2](a) and 12(c)], which 
will be discussed in the next section, the experimental data are found in fair agreement with 
the theoretical predictions. 



C. Specific shear viscosity as a function of T and M 

The BCS pairing gaps and entropies for 88 Mo and 106 Sn are shown in Fig. H]as functions of 
T at several values of M. The general trend of the pairing gap is to decrease with increasing 
T and M. But a slight pairing reentrance occurs for protons in 88 Mo at M = 4 h, and for 
neutrons in 106 Sn at M = 9 and 10 H, where the ga p in creases with T at T < 0.2 MeV. For 
the detail discussion of this phenomenon, see Ref. 13[ and references therein. The pairing 
gap vanishes at M > 22, 18, and 10 h for protons, neutrons in 88 Mo, and neutrons in 106 Sn, 
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FIG. 5: (Color online) Specific shear viscosity rj/s as a function of T at various M (in K) for 88 Mo 
[(a) and (b)] and 106 Sn [(c) and (d)]. The results obtained by using 7/(0) = 1 u are shown in (a) 
and (c), whereas those obtained by using 7/(0) = 0.6 u are shown in (b) and (d) (u = 10 -23 Mev s 
fm~ 3 ). 



respectively at all T. This feature can also be seen in the entropies, which become smoother 
as a function of T with increasing M. Shown in Fig. |5]is the specific shear viscosity rj/s as 
a function of T at several values of M for 88 Mo and 106 Sn. The results are obtained by using 
Eq. (T19]) at two values of r/(0) = lu [Figs. 0(a) and 0(c)] and 0.6m [Figs. 0(b) and 0(d)]. 
The general trend of rj/s to decrease with increasing T at a given value of M is seen in all 
cases. The ratio rj/s also decreases with increasing M at a given T. By using 77(0) = lu it 
is seen that the value of rj/s is always larger than the KSS conjectured lower bound up to 
very high T = 6 MeV and M = 100 and 140 h for 88 Mo and 106 Sn, respectively. However, 
by using the lowest available value of 7/(0) = 0.6tj, extracted from the nuclear fission data, 
one finds the limiting values of M to be around 46 and 55 h, respectively, so that 77/s does 
not violate the lower bound conjecture up to T = 5 MeV, which is the temperature where 
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FIG. 6: (Color online) GDR strength function Sl(uj) from Eq. (|2ip (in arbitrary units) for 88 Mo at 
(T = 3 MeV, M = 41 ft) (a) and (T = 4 MeV, M = 44 ft) (b) predicted by the PDM in comparison 
with the preliminary data from Ref. [1(| . 

the GDR and a nucleus are assumed to cease to exist. These maximal values of M are 
found in reasonable agreement with those obtained in the calculations of the critical angular 



momentum in the entrance reaction channel 



26] and the maximum angular momentum of 



the compound nucleus proposed in Ref. [27 ]. 

Finally, by using the temperature dependence of r]/s and the KSS lower bound conjecture, 
it is possible to examine the reliability of the recent preliminary data for the GDR width in 
88 Mo reported in Ref. jlQj . As a matter of fact, the experimental GDR strength function 
is usually extracted from the GDR photoabsorption cross section, which can be well fitted 
with a Lorentz line shape. The latter is composed of two Breit-Wigner distributions [See 
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f.|a 



Eq. (16) of Ref. [8], e.g.]. In the present case, such Lorentzian-like strength function reads 

S L (u) = -=^—[S(u, E GDR ) - S(u, -E EGDR )} , (21) 

J^GDR 

where S(u, Egdr) is the strength function (fT7j) with the main peak (or the location parame- 
ter of the distribution) located at Egdr, whereas S(u, —Egdr) is obtained from S(uj, E G dr) 
by replacing Egdr with —Egdr- The PDM predictions of the GDR strength function Sl(u) 
in 88 Mo at (T = 3 MeV, Af = 41 ft) and (T = 4 MeV, M = 44 ft) are displayed in Fig. El in 
comparison with the preliminary data of Ref. 10j. The theoretical strength functions have 
been normalized so that the peak of the strength function, obtained at (T = 3 MeV, M = 
41 ft), matches the maximum of the corresponding experimental GDR strength function. 
The figure shows that, while the theoretical and experimental line shapes of the GDR agree 
fairly well at (T = 3 MeV, M = 41 ft) with the FWHM T ~ 11 MeV, they strongly mis- 
match at (T = 4 MeV, M = 44 ft), where the experimental GDR peak becomes noticeably 
narrower with a width T ex ~ 7.5 MeV [See also Figs. [2] (a) and [2] (c)]. By using this value 
T ex and r/(0) = 0.6 u, one ends up with the value of 77/s = 0.85 KSS units. Including the 
error bars in T ex leads to Tf x ~ 6 MeV and T> x ~ 8.5 MeV, which give the values of rj/s 
equal to 0.69 and 0.94 KSS, respectively. All these values are smaller than the KSS lower 
bound conjecture. Therefore, one can conclude that either (i) the data analysis of the GDR 
strength function for 88 Mo at T = 4 MeV and M = 44 ft is inaccurate, or (ii) a violation of 
the KSS conjecture has been experimentally confirmed for the first time ever. A reanalysis 
of the data is expected to clarify which one from these two conclusions holds. 



IV. CONCLUSIONS 



The present paper extends the PDM to include the effect of finite angular momentum on 
the damping of GDR at finite temperature. The formalism is based on the description of 
the non-collective (single-particle) rotation of spherical systems. This implies that the total 
angular momentum J can be aligned along the z-axis, therefore it is completely determined 
by its projection M on this axis alone. 

The numerical calculations were carried out for two spherical nuclei 88 Mo and 106 Sn, for 
which the experimental data for the GDR width are available at T 7^ as well as M ^ 
at sufficiently large values of M. The analysis of the numerical results show that the GDR 
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width increases with M at a given value T for T < 3 MeV. At higher T, the GDR width 
approaches a saturation at M > 60 K for 88 Mo and > 80 h for 106 Sn. However, the region 
of M > 60 goes beyond the maximum value of M, up to which the specific shear viscosity 
rj/s has the values not smaller than the KSS lower bound conjecture for this quantity. This 
maximum value of M is found to be equal to 46 and 55 H for 88 Mo and 106 Sn, respectively, 
if the value r/(0) = 0.6xl0~ 23 Mev s fm~ 3 for the shear viscosity at T = is used. 

A check by using the KSS lower bound conjecture for the specific shear viscosity and the 
same n(0) = 0.6xl0 -23 Mev s fm -3 also shows that the experimental data for the GDR line 

n 

shape in 88 Mo at T = 4 MeV and M = 44 H [lOj leads to a violation of the KSS conjecture. 
This calls for the need of reanalyzing the recent experimental data reported in Ref. 10[ for 
the GDR in 88 Mo at these large values of temperature and angular momentum. 
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